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1. Introduction 

In Quantum Chromodynamics, the current theory of the strong force, many states do not ap- 
pear directly in scattering experiments, but only indirectly in the behaviour of the scattering cross 
sections of observable particles. This is because these states, known as resonances, are unstable 
and decay in a very short time relative to the scattering experiment. Extracting the decay widths 
and masses of these states is thus an important theoretical challenge. 

Lattice Field Theory provides a possible way of extracting these resonance parameters in a 
non-perturbative fashion. The typical method for extracting particle masses in lattice field theory is 
to study the decay of a correlator with the same quantum numbers as the particle in question. The 
large-time behaviour of the correlator is then 

limC(f) =Ze-" lt , (1.1) 

where m is the mass of the lightest particle with the chosen quantum numbers, which can be ex- 
tracted by fitting the correlator. This method however will not work for resonances. By virtue of 
being unstable, resonances are above the multiparticle threshold in their channel and never domi- 
nate the behaviour of the correlator in a simple manner. Fundamentally resonances are not energy 
eigenstates of the Hamiltonian, but rather poles of the S -matrix and so are a truly dynamical phe- 
nomena. 

Given their relation to dynamical scattering processes, resonance parameters can be found 
using the scattering phase shift 8(p). The difficulty lies in obtaining information about 8(p) on 
a Euclidean lattice. In Ref [jl|] it was discovered that there is a connection between the behaviour 
of the two particle energy spectrum in finite volume and the scattering phase shift S(p). Hence 
provided one can accurately determine the energy spectrum it should be possible to obtain 8(p) 
and through it the resonance parameters. Recently another method has been proposed in Ref [0]. 
This method takes the intuition gained from Liischer's method to construct a probability distribution 
which measures the relative frequency of energy levels in the interacting case (resonance present) 
and the noninteracting case (no resonance); we will refer to this method as the histogram method. 
It can be shown that the parameters of this probability distribution are fundamentally related to 
the parameters of the resonance. Furthermore, the method provides a visual tool. The resonance 
should manifest itself as a peak in the distribution. 

Our first discussion on this topic, with an historical introduction and a basic theoretical back- 
ground, can be found in Ref In this work we aim to compare and contrast these two methods. 
In particular we analyse how accurately resonance parameters can be extracted using both methods 
and also the ambiguity in applying the two methods. A first attempt to test the histogram method 
on a simple one dimensional model can be found in Ref [ffl|. 

2. Theoretical background 

2.1 Two particles in a box 

In the continuum, two identical non-interacting bosons of mass m n characterised by a relative 
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momentum p, in a box of volume V = 11;= l ^<> nave a tota l ener gy E given by 



2\lm\ +P 1 



(2.1) 



where due to the finite volume, the momenta p,- are given by pi = jjiii with n,- G Z. On the lattice 
the space-time discretization can have a strong effect (see Figure |l|) in particular when the volume 
is small (large momentum) and m % is big. The correct expression for the simplest discretisation of 
the free scalar field is 

-l 



E = 4sinh 



(2.2) 



where p,- = 2 sin j-rij. It is also valuable to use this expression to describe the energy of interacting 
particles as we will show. 

In a general case such as QCD, where 



K)3.0 




Figure 1: The total energy E for four different levels in the 
continuum (black lines) and in the lattice (red lines) case 
versus L. 



an expression like Eq. ^2| is not available, 5 
we need to determine the non-zero-momentum 
single-particle energy levels numerically and 
then, to determine the two-particle energy 
spectrum, simply multiply the results by a 
factor of two. 

We will focus on the case of a cubic 
lattice, characterised by a single side length 
L; moreover, in a cubic box if n 2 = £? =1 n 2 
is fixed, degenerate energy levels for dif- 
ferent values of n,- can appear. 

In Figure [T] we show a plot of the two 
formulas where it is evident that for small 
volume (L < 15) and a mass m % = 0.46 the 

two spectra are very different; therefore we cannot use the continuum formula to describe our 
Monte Carlo results. 

2.2 Avoided level crossing 

Let us introduce another particle a in the box (at the moment, not interacting) with mass m a ; 
we are interested in studying the elastic scattering between the n particles therefore we impose 
the constraint 2m n < m a < Am n . In Figure § (Left) the a energy level is the horizontal line that 
intersects the two-particle levels at various system sizes L. 

In Minkowski space if we introduce a three point interaction okk between the fields, the a 
can be unstable and decay into two % -particles. In Euclidean space and in a finite volume the 
scenario is different; because of the interaction, the energy eigenstates are a mixture of a and 2% 
Fock-states and they are all stable. The mixing is manifested in avoided level crossings (ALCs) in 
the energy levels as shown in Figure ^ (Right). 

2.3 Luscher's method 

The best known method for analysing resonances was proposed by Luscher (Ref [Q]). This 
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Figure 2: (Left) The spectrum of a system of two non-interacting particles of mass m K = 0.4544 worked 
out using Eq. 2.2; the horizontal line describes the particle a at rest with mass m a = 1.3517. With these 
parameters the intersection between a and the two-particle level n 2 — 1, i.e. (1,0,0), is set at L = 12. (Right) 
Avoided level crossings where on the (Left) there were intersections between a and 2%. 



method involves using information on the scattering phase shift. Since the scattering phase shift 
contains information on resonance parameters this provides a way to extract the resonance mass and 
width. The scattering phase shift itself can be obtained using the relationship found in (Ref [[I]]). 
In essence, this relationship provides a mapping between the values of the two-particle energy 
spectrum in finite volume and the scattering phase shift in infinite volume. 

The relationship is proven first in non-relativistic quantum mechanics It then holds in quantum 
field theory as the relativistic case can be cast in a non-relativistic form, with the Bethe-Salpeter 
Kernel playing the role of the potential. This is achieved by means of an effective Schodinger 
equation, first constructed in (Ref |^]). The precise relationship is 1 



tan(>(K)) 




<$>{K) + Kn , 



K ■ 



pL 

2k ' 



(2.3) 
(2.4) 



where p is the relative momentum between the two pions. ^j s (l;q 2 ) is a generalised Zeta function, 
given by 



where Yj s {d,§) are the spherical harmonics. Eq. |23]is known as Liischer's formula. 
To obtain resonance parameters using this relationship one proceeds as follows: 



(2.5) 



1. Through Monte Carlo simulations, obtain the two particle energy spectrum E n (L) at different 
volumes; 



2. Through dispersion relations obtain a momentum from the energy spectrum, p n (L); 



'in truth the relationship is more general than this, the formula quoted here is for the spin-0 channel, which is the 
only one relevant here. 
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3. Through Eq. 2.3 p n (L) gives a value for 8(p n (L)); 

4. By repeating this procedure for several energy levels and volumes, one obtains a profile of 

8(p); 

5. This profile of 8(p) is then fitted against the Breit-Wigner form for the scattering phase shift 
in the vicinity of a resonance: 

This fit should give the resonance mass M a and width F a . 

The work in this paper applies this method to the the 0(4) model in the broken phase and also 
compares its performance against a more recent proposal. 

2.4 Histogram method 

An alternative method to determine the parameters of a resonance is based on a different way 
to analyze the finite volume energy spectrum (Ref [§]). The basic idea is to construct the probability 
distribution W(E) according to the prescriptions: 

1. Measure the two-particle spectrum E„ (L) for different values of L and for n = 1 , • • • ,N levels; 

2. Interpolate the data with fixed n in order to have a continuum function E n (L) in an entire 
range Le [L ,L M ]; 

3. Slice the interval [Lq,Lm] into equal parts with length AL = (Lm — Lq) /M; 

4. Determine E n (Lj) for each L, (i = 0, • • • ,M); 

5. Choose a suitable energy interval [E m i n ,E max ] and introduce an equal-size energy bin with 
length AE; 

6. Count how many eigenvalues E n (Li) are contained in each bin; 

7. Normalize this distribution in the interval [E m i n ,E max \. 

Actually, the distribution considered in Ref [^] is W(p) but this does not have an important 
effect on our analysis; as a matter of fact, the relation between them is: 

W(p)=W(E)(^\ , (2.7) 



where the correct dispersion relation we should use is Eq. |2.2fc the multiplicative term will not 
modify the Breit-Wigner shape near the resonance. 

It is possible to show that the probability distribution W(p) is given by W(p) = cj^ = i \p' n (L)]~ 
and differentiating the Luscher formula with respect to L, it turns out (c is a normalization constant): 



n P ) = - 1 

V n=\ 



c ^ T . . 2x8' {p) 
Ln(p) + - Fl 



§'{<ln{p)) 



(2.8) 
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The authors of Ref [g] introduce Wq(p), which is determined by Eq. 2.8 with 8{p) = and L n {p) 
corresponding to the free energy levels: they show that in order to subtract the background (free % 
particles) it is convenient to consider the subtracted probability distribution W(p) = W{p) —Wo(p). 

Using convenient approximations in the Luscher formula and in the limit of infinite number of 
energy levels (infinite volume) it turns out: 

W( P ) - W (p) « - ( ^ - 8'(p)) . (2.9) 

p V p J 

This last quantity is determined by 8(p) alone and close to the resonance, assuming a smooth 
dependence on p for the other quantities, it follows the Breit-Wigner shape of the scattering cross 
section with the same parameters: 

>nri-w >(p)2 _4 |2+Mgr2 . (2.io) 

In Ref [||] this method is tested on synthetic data produced using the Luscher formula by 
experimentally measured phase shifts. The main task of our work is to test this method on an 
effective field theory where a resonance emerges, producing data by lattice simulations. 



3. The model 

The model we have used in our simulations is essentially the 0(4) model in the broken phase. 
This model has previously been used to test Liischer's method (Ref |Q]). The Lagrangian is the 
following: 

£? = ~d<j>id<l>i + l(<l>?-V 2 ) 2 -mlv<!>4 , with i=l,2,3,4 . (3.1) 

The term proportional to 04 is introduced to break explicitly the symmetry and therefore to give 
mass to the three Goldstone bosons. To understand the meaning of the terms and the parameters in 
the Lagrangian, we first introduce the new fields a and p,- (with the constraint p;p, = 1): 

<t>i = {v + a)pi, with i = 1,2,3,4 ; (3.2) 

then, we expand the potential around the classical minimum = v 2 (using also p,dp,- = 0): 

1 ? 1 ? 1 

££ = -v z dpidpi + -(J z dpidpi + -dada + vadpidpi 

+ Aa 4 + 4v 2 Aa 2 +4vAa 3 - m 2 v 2 p 4 - m 2 vap 4 . (3.3) 

The field a is clearly related to the massive field whereas the four constrained fields p, are related 
to the three "pions". There is an easy way to see this based on the treatment of the non-linear sigma 
model 2 ; we introduce the pions using an element of SU (2): U = exp (jjCjGjj, where Oj are the 

2 See for example Ref |||| Sec 2.4.3. 
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three Pauli matrices and / is the pion decay constant. It turns out that 

z 1 j=i 



(3-4) 



On the other hand, we can introduce an element of SU (2) by U = p 4 + iOjPj, with j = 1,2,3 and 
the constraint p| + pjpj = 1; in this case we have: 



1 3 
-Tr (dpUd^W) = dp 4 dp 4 + £ dpjdpj = £ apidp,- 

Z 7=1 



4 

I 

1=1 



Comparing Eq. 3.4 and Eq. 3.5 it turns out (f — > oo): 



(3.5) 



4 1 3 

£ dp/dp,- - £ d^dn/ • ( 3 - 6 ) 
i=i / j=\ 

Using the same previous approach and the relation Tr(£/ + £/') it is easy to show that 

p 4 = -^7ij7ij + const (3.7) 
We can rewrite the Lagrangian of Eq. 33 using Eq.s ^^ - |3.7| and introducing ftj = Ujj : 



1 m 2 
+ Xo A + Av 2 Xo 2 + AvXo^ + -m 2 n Ttjftj + ^oftjftj . 



(3.8) 



In this expression, it is now evident that Ttj are the pions with mass m n and a is a massive field 
with mass m a = 2vy2X. It is interesting to note there are two three-point interaction terms, both 
inversely proportional to V. 



4. Monte Carlo simulation 

The theory described by the Lagrangian Eq. |3j] was simulated using an overrelaxation algo- 
rithm for the first three fields followed by a Metropolis update to guarantee the ergodicity and a 
Metropolis algorithm for the field ty 4 . 

4.1 Single and two-particle spectrum 

In order to determine the single particle spectrum we first introduce the partial Fourier trans- 
form (PFT) of the four fields 0,-: 

1 — 2k 

= y L & t) e ~ aP > Pi = -T n i > n i = °> • • • > L * - 1 • ( 4 - 
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The single particle mass is extracted from the zero momentum correlation function (n = 0): 

Q(t) = (4>i(n,t)U-n,0)) , (4.2) 

in particular with i = 1 , 2, 3 we can determine the mass of the %i particles; with i = 4 we extract the 
mass of the a particle. Because of the different way we update the four fields, it turns out that the 
mass m % is determined with a higher precision then m a ; actually, this is not a real problem because 
we are mainly interested to a good determination of the "pion" mass. 

The two-particle spectrum is measured by introducing operators with zero total momentum 
and zero isospin: 

3 

o*(0 = E&(n.0&(-ff>0; (4-3) 

i=l 

we take in account N—l different operators, corresponding to n 2 = 0,\, - ■ ■ ,N — 1 . A N-th operator, 
that clearly has the correct quantum number, is the PFT of the field a (actually $4) with p = 0. 
To determine the energy levels we use a method, introduced in Ref [B], based on a generalized 
eigenvalue problem applied to the correlation matrix function Qj(t) = (0,-Oj), i.e. a matrix whose 
elements are all possible correlators between the N operators. 

4.2 Numerical results 

In order to test the applicability of the two methods for different widths of resonance, we 
consider three different sets of parameters. The first set is characterised by v = 1.0, X = 1.4, 
m K = 0.36. We tuned these parameters to have the intersection between the a energy level and 
(1,0,0) two-particle energy level around L = 12. The physical mass for the pion turns out to be 
m Tc h = 0.460(2). The spectrum (the first 6 levels) was determined for different volumes (8 < L < 19) 
and the relative error works out to be in the range 0.5% - 1.0% (see Figure || (Left)). 
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Figure 3: (Left) Spectrum of the theory for different values of the volume for the following simulation 
parameters: V = 1.0, A = 1.4, m K = 0.36. The dashed lines describe the free two-particle spectrum. (Right) 
The interpolated data using a polynomial. 

First of all, we interpolate the data for each level using three different polynomials of order 3, 
4 and 5 in order to have a relation between the energy E and the side box L in the entire interval 
[8, 19] and to provide a way to evaluate the systematic errors in our final results; in Figure |3| (Right) 
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we show one of these interpolations. In Figures [| the dashed lines describing the free two-particle 
spectrum are calculated using Eq. Here we note that the mass m n in Eq. 22 is not the bare 
mass. The value measured on the lattice m^ 1 = 0.460(2) is taken as the physical value; the relation 
between them is well-approximated by m^ 1 = 2sinh~ 1 (m 7r /2). 
This free spectrum is used to deter- 
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mine the distribution Wq which is then sub- 
tracted to W that is obtained from the in- 
teracting spectrum. It is important to note 
that if the number of levels considered to 
plot W in the interacting spectrum are N, 
then the number of levels we have to con- 
sider in the free spectrum to determine Wo 
are just N — 1. 

Using the three previous polynomials 
we are able to produce a large number of 
data (we fix AL = 0.001) that we can then 
use to get the probability distribution W 
described in Sec |Z4] with the correspond- 
ing systematic errors; fixing the bin width to AE = 0.005 we obtain the histogram W of Figure |j. 
Note that to get W both W and Wo are worked out from the same range with L G [8, 19]. The er- 
ror bars in Figure |] are the results of the systematic errors coming from the histogram W and the 
statistical errors coming from the histogram Wo- 

Clearly, the shape of the histogram in Figure ^is far from the Breit-Wigner shape; the reason is 



1.4 
mass 



Figure 4: The probability distribution W — W — Wq ob- 
tained by data from Figure [3| 



related to the fact we are considering only six energy levels but the conclusions of Sec 12^ are true 
only in the limit of an infinite number of levels. Moreover a lot on jumps and spikes are present. 
Our task is now to try to improve this result in order to get more information from our raw data. 

We investigated the origin of the spikes 
and we understood they are related to a 
"wrong" background Wo. It is easy to see 
that the spikes appear every time there is 
the intersection between the six levels of 
the interacting spectrum (or of the five lev- 
els of the free spectrum) with the extrem- 
ities of the volume range (L = 8 and L = 
19). Near those two extremities we have 
to be careful with what is the correct back- 
ground; Figure || shows a corrected back- 
ground subtraction. In order to correctly 
subtract the free background, we lengthen 
each free spectrum line. This is done so 

that the extremity of that line has an energy equal to that of the extremity of the interacting spectrum 
line closest to it. In this way all interacting lines are subtracte correctly rather than the subtraction 
being affected by the limit of the volume range that we are actually using in our simulations. Using 




Figure 5: Energy levels of Figure |] (Right) with the cor- 
rect free two-particle spectrum background. 
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this procedure to determine Wo we get the correct histogram of Figure ^ (Left). Unfortunately, in 
Figure ^ (Left) we continue to see a jump for E rj 1.35; the origin of it can be understood looking 
at Figure ||. There are two extremity lines, one at L = 8 and one at L = 19 (both around E rj 1.35), 
that are without a "background"; actually, in this case the background is the resonance itself we are 
looking for. Therefore, there is no way to avoid the presence of this jump because we do not know 
anything about the resonance; the only thing we can do is to completely exclude from our analysis 
those two levels, hoping that the resonance can appear. In Figure ^ (Right) we show the probability 
distribution W in this last case; now clearly a Breit-Wigner shape appears. 

It is now possible to fit these data to Eq. 2.10| to determine the parameters of the resonance; 
applying a sliding window procedure around the peak, they turn out to be: M a = 1.330(5) and 
r CT = 0.10(5). We have simulated the theory with a second set of parameters, corresponding to a 




mass 



le+02 




1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 



Figure 6: (Left) Probability distribution W obtained by data from Figure pi (Right) Probability distribu- 
tion W obtained excluding from the analysis the two levels that in Figure g are without a corresponding 
background. 

larger width: v = 1.0, A = 4.0, m % = 0.56. In this case, we tuned them to have the intersection 
between the a energy level and (1,0,0) two-particle energy level around L = 8. The physical 
mass for the pion turns out to be m 1 ^ 1 = 0.657(3). In Figure ^ (Left) we plot the spectrum for 
6 < L < 20 for the first six levels; the relative error varies in the range 0.05% - 0.2%. If we 
repeat all the procedure as described before (in particular we exclude the two levels which are 
"without" background) we get the histogram of Figure ^ (Right); also in this case we can clearly see 
a Breit-Wigner shape and we can fit these data obtaining the following parameters: M a = 2.01(2), 
r ff = 0.35(10). 

Finally, we have run a third series of simulations with parameters v = 1.0, X = 200.0, m n = 
0.86. They have been tuned to have the intersection between the a energy level and (2,0,0) two- 
particle energy level around L = 10. Because in this case we are considering higher momentum, 
we expect the width of the resonance is larger then the previous cases. In this case we take in 
account 13 levels to describe better the shape of the resonance. In Figure || (Left) the spectrum 
for 6 < L < 15 is plotted; the relative error varies in the range 0.15% - 0.4%. The physical mass 
for the pion turns out to be m^ h = 0.938(3). Unfortunately, as it is shown in Figure [8] (Right) the 
probability distribution plot is flat, i.e. no Breit-Wigner shape emerges. It is clear that in this case, 
the only way to determine the parameters of the resonance is to increase considerably the number 
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1.7 1.8 1.9 



2.1 2.2 
mass 



2.3 2.4 2.5 2.6 



Figure 7: (Left) Spectrum of the theory for different values of the volume for the following simulation 
parameters: v = 1.0, A = 4.0, m K — 0.56. (Right) Probability distribution W using the correct background 
and excluding the two levels that are without a corresponding background. 
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Figure 8: Like Figure [7| but with simulation parameters: v = 1.0, A = 200.0, m K = 0.86. 



of measurements and consequently to decrease the relative errors in the spectrum determination. 
In Table [T] a summary of our results for the three sets of parameters are shown. 



Relative error in E(L) 




8{M a )/M a 


r ff 


8{T a )/T a 


0.5%-1.0% 


1.330(5) 


0.4% 


0.10(5) 


50% 


0.05%-0.2% 


2.01(2) 


1.0% 


0.35(10) 


28% 


0.15%-0.4% 











Table 1: Results for the three sets of simulation parameters with the corresponding relative errors. 

In applying Liischer's method to the data one only needs the original raw data; there is no need 
to fit it to a polynomial expression as in the histogram case. The first step is to convert the data 
on the energy spectrum to data on the momentum spectrum. This requires the dispersion relations. 
However should it be the lattice or continuum dispersion relations? The lattice dispersion relations 
are more natural, since they suppress lattice artifacts, but results were obtained for both below to 
emphasise how much more effective they are. In order to use p n (L) and Eq. to obtain 8(p), 
knowledge of 0(fc) is needed. The main difficulty here is the cumbersome definition of £¥j s (\;q 2 ). 
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Figure 9: (Left) 8(p) using Lattice dispersion relations at: v = 1.0, A = 1.4, m K = 0.36. (Right) Same 
parameters, but with continuum dispersion relations. Both done with our approximation. 



Not only that, but as mentioned the summation expansion given above does not even converge in 
the region required. A more convenient definition is an integral representation of 3?j s (l;q 2 ) given 
in Appendix C of Ref |jl]]. This expression can be used to numerically evaluate ^j s (\;q 2 ). Some 
data on the values of (?c) can then be obtained from Eq. |2~4| . We fitted (j) ( k) against these values 
to obtain an approximation of 



» (-0.09937) K 8 + (0.47809) k 6 
+(-0.62064) k 4 + (3.38974) k 2 

The error between this approximation and 0(jc) is negligible compared with other errors. 



(4.4) 



From here we use Eq. 2.3 to obtain a profile of 8(p). For the narrow case one can see the 
difference between use of the continuum dispersion relations and the lattice dispersion relations 
in Figure |9[ The lattice dispersion relations provide a tighter fit of the data, as well as having 



Results 


Parameters 


m 




v = 1.0, A = 1.4 


M a = 1.35(2) 


M ff = 1.36(4) 




T ff = 0.115(8) 


r ff = 0.17(2) 


v = 1.0, A =4 


M a = 2.03(2) 
F a = 0.35(2) 


M a = 2.2(2) 
r ff = 0.42(5) 


v = 1.0, A =200 


M a = 3.1(7) 


M a = 3(1) 




r a = 1.2(5) 


r ff = 2(i) 



Table 2: Resonance mass and decay width using two different approximations for <P(k). 

smaller errors. We also compared the use of the traditional approximation of 0(k) = kk 2 with 
our approximation. After fitting, the results for the resonance mass and decay width in the two 
approximations are (both using lattice dispersion relations) shown in Table ||. 



The errors are smaller when the approximation of Eq. 4.5 are used, particularly for the broad 
resonance. It should also be noted that the two approximations effect the two resonance parameters 
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differently. The dispersion relations have a more direct effect on M a while the approximation of 
<P(k) has a greater effect on r a . This is because a cruder approximation of 0(k) effects the profile 



Results 


Parameters 


Luscher's Method 


histogram method 


v = 1.0, A = 1.4 


M a = 1.35(2) 


M c = 1.33(5) 




T CT =0.115(8) 


r a = 0.10(5) 


v = 1.0, A = 4 


M a = 2.03(2) 


M ff = 2.01(2) 




r ff = 0.35(2) 


T a = 0.35(10) 


v = 1.0, A =200 


M ff = 3.1(7) 


M ff =N/A 




r ff = 1.2(5) 


r ff = w/A 



Table 3: A comparison between the Liischer and the histogram method. 



of the scattering phase shift, which is related to the decay width. These results suggested it is 
optimal to use the lattice dispersion relations and our approximation. 

4.3 Comparison between the two methods 

The results for Luscher's method compared with the histogram method are shown in Table |3| 

Luscher's method gives smaller errors 
than the histogram method, but the results 
are broadly consistent. Luscher's method 
manages to provide some estimate on the 
width of the resonance in the broad case. 

The broad resonance becomes a prob- 
lem for the histogram method because there 
is no obvious peak to indicate the reso- 
nance mass (and hence no width of that 
peak to determine the decay width). One 
would need very precise data in order to 
avoid a washing out of the structure of the 
histogram. Luscher's method also becomes 
more difficult to apply in the case of broad 
resonances. In the case of a broad reso- 
nance the profile of 8{p) is quite flat, hence a large range of parameters will be capable of fitting to 
the profile. Again an accurate determination of the energy levels is required to determine the pro- 
file precisely enough so that this is prevented. Considering the amount of work necessary before 
one can use the histogram method (as detailed above), Luscher's method is considerably easier to 
apply, provided one has a good approximation of <j>(K). However, the histogram method can be 
used as a visual tool for spotting the resonance. 

One restriction of Luscher's formula is that it only applies in the elastic region. An example 




Figure 10: Inelastic data with Luscher's formula. For 
the case of v = 1 .0, A = 1 .4. (Onset of inelastic region 
marked). 



of what happens in the inelastic region is provided in Figure 10. It is possible that the histogram 
method will provide a means of determining the presence of a resonance in the inelastic region. 
Certainly a histogram can be constructed in the inelastic region, the only difficulty is that with the 
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inapplicability of Luscher's formula it is unclear that the parameters of this histogram will have 
any relation to those of the resonance. 

5. Conclusions 

We have compared and contrasted the Luscher and the histogram methods. Luscher's method 
appears to both easier to apply and give smaller errors, however the histogram method does give 
results consistent with Luscher's method and does indeed visually indicate the presence of a reso- 
nance. 

There are two major difficulties with both methods. First for broad resonances the relevant 
structure is washed out to some degree. For a histogram, the peak is hard to locate, while for the fit, 
the profile of the phase shift is poorly constrained. Secondly there is the inelastic region. Luscher's 
formula cannot be used there. The histogram method can be applied to the data, but there is no 
argument that this is a sensible thing to do. There is also a difficulty in the general case, relevant 
to QCD, which has not been examined here. In the model above the resonance is clearly present in 
the channel, since this is an explicit feature of the Lagrangian of the model. In general however a 
resonance may not be so obvious and there is no reason a priori to expect that it will have a purely 
Breit-Wigner form. 
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